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I. INTRODUCTION 

The past two decades have produced impressive 
progress in the reliable preparation and manipulation of 
single quantum systems. To achieve this, one of the ma- 
jor challenges has always been the isolation of individ- 
ual systems and the suppression of decoherence. Over 
the years, the efforts in this direction moved from simply 
perfecting technical aspects of experiments that minimise 
environmental effects, to more elaborated strategies to 
preserve quantum coherence. Several tactics have been 
considered in this context, including reservoir engineer- 
ing [U [2] , dynamical decoupling [3] H] , decoherence- free 
subspaces (DFS) [5 , 6 , open loop control [7 and feedback 
methods [MO]. 

Many strategies to control decoherence have been de- 
veloped recently in the context of quantum information, 
with a focus on the preparation and preservation of en- 
tangled states. Schemes to protect entanglement using 
reservoir engineering [11] , quantum feedback [T2HT8] , and 
reservoir monitoring [19] have been proposed, while long- 
lived entangled states have been produced in the lab us- 
ing DFS ideas [20] [21]. A remaining challenge in this 
field is to devise methods to control entanglement in large 
quantum systems. As the number of components grows, 
in general the fragility of entanglement against decoher- 
ence increases [22H24] and the design of control methods 
becomes more demanding. 

In this paper we address the problem of stabilising 
multipartite entanglement using direct, or Markovian, 
feedback based on continuous monitoring of quantum 
jumps. In this case, feedback is directly proportional 
to the measurement signal and a master equation repre- 
senting the unconditioned dynamics of the system can be 
obtained [8]. We show, using an extension of the model 
proposed in [16] [17] for multiple qubits, how the intro- 
duction of feedback can change the stability of the steady 
states of the uncontrolled master equation and be used 
to protect highly entangled states. In this case, feedback 
is used to engineer the reservoir by modifying the deco- 
herence operator, which selects preferred steady states 



of the system. In particular, we show a general strategy 
to generate and stabilise specific entangled steady states 
and describe a simple and realistic feedback control for 
the case of four qubits. 

The paper is organised as follows. In Section [TT| we 
briefly discuss the general structure of the master equa- 
tion describing a system under jump-based feedback and 
its consequences to the properties of steady states. In 
Section |III| we describe the model where entanglement 
generation will be investigated. In particular we analyse 
how the steady state entanglement in the system scales 
with the number of particles. The effects of feedback in 
the structure and dynamics of subspaces is presented in 
Section IV In Section[V|we describe a general strategy to 
manipulate the steady states of the system to protect spe- 



cific target states. In Section [VT] we apply this strategy in 
the case of four particles and present an alternative sim- 
ple local feedback to produce stationary entanglement. 
Section [VII] finishes with a discussion about the results 
and future perspectives. 



II. STEADY-STEATE PROPERTIES OF 
QUANTUM-JUMP-BASED FEEDBACK 

The evolution of an open quantum system interact- 
ing with a single-channel environment can be described, 
under the assumption of complete positivity and Marko- 
vian dynamics, by a master equation in the Lindblad 
form mi [261: 



P = l v \c\p = 7 



cpc 



(1) 



where 7 measures the strength of the system-environment 
coupling and c is the operator that describes the effect of 
the bath on the system. 

To discuss feedback dynamics, we first need to consider 
how the system is being measured. Eq. ([TJ corresponds 
to the dynamics of an ensemble of quantum systems de- 
scribed by the density matrix p and does not convey 
any information on how a single quantum system would 
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evolve under measurement. Consider instead the situa- 
tion where information about the system is extracted by 
continuously monitoring the environment, more specifi- 
cally in terms of quantum jumps. The time evolution 
is then conditioned on the measurement outcomes, and 
can be described in terms of stochastic quantum trajec- 
tories [271 HE]- The dynamics will intersperse random 
quantum jumps, defined by the action of the operators c, 
with periods of continuous evolution, given by the non- 
hermitian effective Hamiltonian H e ^ = —ifvyc^c/2, cor- 
responding to periods with no detection. 

Now consider that the detection events trigger feed- 
back pulses, which are assumed to be short compared 
to the other timescales in the system. In this case, 
the action of the jump operator c is followed by the ef- 
fect of feedback and the state transforms according to 
— » Uc\ip), where U is a unitary operator representing 
the feedback [8 . In the absence of detection, no feed- 
back is applied and the system evolves with the previous 
effective Hamiltonian for the no-jump evolution. 

Because the feedback is assumed to be instantaneous, 
one can also obtain a Markovian master equation for the 
jump feedback [8 . Following the reasoning above for a 
single trajectory, one can see that the addition of feed- 
back corresponds to simply replacing the jump c by the 
combined action Uc in the dynamics. The master equa- 
tion then becomes [8] 



p = jV[Uc\p. 



(2) 



Feedback can then be viewed as an instance of reservoir 
engineering as it allows us to modify the decoherence 
operator in the master equation through the unitary U. 
This modified irreversible evolution can then be used to 
dynamically protect quantum states. The idea is to en- 
gineer the steady states of Eq. Q such that the system 
naturally evolves towards the desired target state. We 
therefore need to analyse the stationary states of Eq. ^ 
and their stability conditions. The mathematical aspects 
of the stability of steady states of Lindblad master equa- 
tions have been addressed in the past [29-31, including 
the application in Markovian feedback using homodyne 
detection [32]. Here we will only briefly discuss some of 
the basic properties that will allow us to construct feed- 
back strategies based on quantum-jump monitoring. 

Note first that the pure steady states, \ip 3 s)i °f Eq. ([2| 
are the eigenstates of c with eigenvalue [43]. Note also 
that the introduction of feedback does not affect these 
states since we also have Uc\i/} ss ) = 0. The dark states 
are therefore uniquely determined by the structure of the 
jump operator c and not by the feedback. Although 
this limits our ability to engineer the dynamics using 
quantum-jump control, feedback can still be extremely 
useful to modify the stability properties of the steady 
states. Indeed, even though the combined jump Uc can't 
move the system once it is in a steady state, it can in- 
duce transitions from outside into the dark subspace and 
be used to select a single dark state to the exclusion of 
others. 



In the remainder of this paper we will analyse this gen- 
eral jump- feedback structure in the context of the model 
proposed in [16] with the goal of generating multipartite 
entangled steady states. 



III. MODEL 

The system consists of TV two-level atoms resonantly 
driven by a laser with Rabi frequency ^, and coupled, 
with strength g, to a single cavity mode as depicted in 
Fig.[l] We assume that the light field leaks through one 
of the cavity mirrors at a rate k and is detected by a 
photodetector. When the photodetector registers a pho- 
ton, a short feedback pulse is applied to the atoms. In 
the limit where the cavity decay rate is large compared 
to the other rates of the system, we can adiabatically 
eliminate the cavity mode [TSJ [17] and obtain the master 
equation: 



-in 



(j+ + 



(3) 

where T = g 2 / n is the effective collective decay rate. 
The angular momentum operators are defined as J± = 
a j,±i where aj- = \gj)(ej\ are the lowering operators 
for the atomic electronic levels. We also consider the 
possibility of atomic spontaneous emission with rates jj , 
as described by the third term in the master equation. 
At this stage we will not assume any specific form for the 
feedback Hamiltonian i^fb, only that its net effect is the 
action of the unitary evolution 



U = exp iH^St/h 



(4) 



The second term in Eq.([3| is of the form of the feed- 
back master equation described in the previous section. 
As discussed before, the steady states from this part of 
the dynamics will be given by the eigenstates of the oper- 
ator J_ with eigenvalue zero. However, these states won't 
be steady states of the full equation since we have extra 
Hamiltonian and decay terms. Ignoring the spontaneous 
emission effects for the moment, we can see that simulta- 
neous eigenstates of J_ and J + with eigenvalue zero will 
not evolve under the Hamiltonian nor under the feed- 
back terms, and will determine the subspace of steady 
states of the system. In what follows we will describe 
this subspace structure and the entanglement properties 
of steady states. 



A. Subspace structure 

The Hilbert space of a two-level atom is equivalent to 
that of a spin- 1/2 system. Our system then consists of 
N of these pseudo-spins and it is convenient to use an 
angular momentum basis to describe the problem. The 
eigenstates of this multipartite system can be found by 
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FIG. 1: Diagram of the model. A series of N two-level atoms 
are coupled resonanltly to a cavity and equally driven by an 
external laser field. Light that leaks from the cavity is de- 
tected in a photodetector, triggering a feedback pulse. 



the usual rules of addition of angular momentum and the 
relevant quantum numbers are the total spin J and the 
total spin projection J z . The dynamics in our system is 
dictated by the operators J + and J_, whose action on 
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FIG. 2: N- concurrence for various kinds of states as a function 
of the number of atoms. Vertical bars indicate the range of the 
concurrence for the many-body singlets forming the subspace 
of pure steady states of Eq. (J3| with 7* = 0. Concurrence 
for GHZ, W and linear cluster states are shown by dashed, 
solid, and dotted lines, respectively. Solid squares represent 
the stationary states of Eq. ^ in the absence of feedback, 
when all atoms start in the ground state. 



the eigenstates \J,J Z ) is given by 



J±\J, Jz, Aj) = V(JT J z ){J±J z + l)\J, J z ± 1, Aj). 

(5) 

Within each total spin subspace (excepting where J = 
N/2) there exist smaller subspaces that are closed under 
J±, and these are distinguished by the extra index Aj. 
The index J z determines different states within these sub- 
spaces that can be connected through the J± operators. 
The division of the full J subspaces into these smaller 
subspaces is not unique, but our analysis is independent 
of the choice of basis that defines this division. Note also 
that the subspace of pure steady states discussed before 
corresponds to J = 0, a situation only achievable for an 
even number of particles where J can assume the values 
0,1,...,^. These states, known as many-body singlets, 
have interesting entanglement properties [33 - 36 with ap- 
plications to quantum communication [37] and computa- 
tion [38] . 



B. Multipartite entangled steady states 

Having described the subspace of pure steady states, 
we are now in position to investigate the entanglement 
properties of all dark states within this subspace. To do 
this we will use the Cn concurrence [22 , 39 as a measure 
of entanglement in the system. Although the correlations 
in a TV-partite system cannot be completely characterised 
by a single scalar quantity, the Cn concurrence is useful 
to analyse the scaling properties of entanglement [22] . 

Let us start with the example of a 4-partite system. In 
this case the subspace of dark states can be parametrised 



as 

|VW = -^=[ae i9 \ggee} + l3e^\gege} 

_ (ae ie + fie*) \geeg) - (ae ie + fie*) \egge) 
+ f3e i *\egeg)+ae i0 \eegg)}, (6) 

where N is a normalisation factor. The TV-concurrence 
of elements of this set varies from \/7/2 to V% depend- 
ing on the values of a, /?, and <\>. This range is plotted 
in Fig. [2] together with the equivalent range for N = 6 
and N = 8 atoms. Curves for representative states like 
GHZ, W, and linear-cluster states [42 are also displayed 
for comparison. These results show that, in the absence 
of spontaneous emission, the model described by Eq. ([3| 
has highly entangled steady states , but it does not auto- 
matically follow that any initial state will lead to a highly 
entangled steady state. In the absence of feedback, each 
component of the initial state with a well-defined angu- 
lar momentum stays within that subspace. If the system 
starts in a dark state then it will remain in that state 
indefinitely. However, if the initial state has a little over- 
lap with the pure dark states, then the system will likely 
evolve towards a mixed steady state with low entangle- 
ment. This situation is represented by the filled squares 
in Fig. [2] which correspond to the steady state concur- 
rence achieved after evolving the system from all atoms 
starting in the ground state. This shows that despite the 
existence of pure steady states, these are not attractive 
points of the dynamics and therefore the system does not 
evolve naturally towards them. Our goal here is to design 
a feedback that could change this scenario by making the 
system dynamically evolve to a target highly entangled 
state. 
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IV. DYNAMICS BETWEEN SUBSPACES IN 
TERMS OF QUANTUM TRAJECTORIES 

The feedback design relies on understanding the de- 
tailed dynamics of the system and its consequences to the 
structure and stability of steady states. A natural way to 
approach this problem is to use the quantum trajectory 
method described briefly in Section [TTJ The method not 
only provides a faithful picture of what a single run of an 
experiment would look like, but also provides a sensible 
splitting of the dynamics in terms of jumps, when con- 
trol acts, and a continuous evolution independent of the 
feedback. 



A. The action of jumps 

When a detection event occurs, the state of the sys- 
tem changes with the application of the combined jump 
operator U J_ : 



B. Interplay between jump and no- jump evolutions 



UJ. 



(J+J-) 



(7) 



The J_ operator does not mix subspaces with different 
J's but rather induces transitions towards the states with 
minimum spin projection (J z = —J) within each sub- 
space. However, the subsequent application of the feed- 
back unitary can cause transitions between different sub- 
spaces and be used to help move the system to a target 
state. This can happen, for example, if the intermediate 
state |^), which results from the action of J_ on the ini- 
tial state | ^0)5 is connected to the targeted steady state 
|^t) via the feedback operator J7, i.e. , 



\ 1 p t ) = u\ 1 p i ) = uj_\ 1 p ), 



(8) 



where normalisation factors have been omitted for sim- 
plicity. Once the system is fully in the steady state then 
no more jumps will happen since UJ-\ip t ) = 0. 

Note, however, that the jump dynamics is not enough 
to account for a full transfer of population to one of the 
pure steady states of the system. First, because states 
with J z = — J are also dark to the action of U J_ and any 
initial population in those states cannot be directly con- 
nected only through jump dynamics to the target states. 
Second, because if the system is not entirely in the dark 
subspace, then a jump event will have exactly the oppo- 
site of the desired effect. Consider, for example, the sys- 
tem initially in a mixed state p = \a\ 2 \<j))((j)\-\-\P\ 2 \tlJt}(^t\^ 
where \<jj) is not in the subspace of steady states. Then, 
the action of the operator U J_ will lead to the final pure 
state = UJ-\<j)) and, unless \<jj) is one of the states 
connected to the target state via the jump operation, 
then the system will be completely removed from the 
dark subspace. This process, however, has a probability 
to occur that goes to zero as /3 — >• 1. 



The dynamics between detection events is given by 
d , 



dr 



+ j_)|^)-ri + j_ 



(9) 
(10) 



where the first term represents the driving Hamiltonian 
of Eq. and the second term is the non-Hermitian ef- 
fective evolution described in Section HI The solution at 
any given time can be found by integrating this equa- 
tion and then renormalising the state. This last step, 
necessary due to the norm decay induced by the effec- 
tive Hamiltonian, is crucial to understand the transfer of 
population in the absence of detections. 

The norm of different components of the state will de- 
cay at different rates, and renormalisation will transfer 
population from states that decay quickly to states that 
decay slowly. Physically, this can be understood from the 
dynamics of the system conditioned on the measurement 
results. As we continuously monitor the system and no 
jump events are observed, we acquire information that 
it is more probable that the system is in a state where 
jumps are less likely to occur. To make this statement 
more quantitative, let us look at the action of the effec- 
tive Hamiltonian on our basis states: 

J+J.\J, J z ,\ k j) = (J + J Z )(J-J Z + l)\J,J z ,Xj) (11) 

The norm decay will be proportional to the factor (J + 
J Z ){J — J z + 1), which is zero for states with J z = — J 
and is at least equal to 2 J for all other eigenstates. This 
means that the renormalisation step will gradually trans- 
fer the population not only to the dark subspace with 
J = 0, but to all the states with the lowest spin projec- 
tion [J z = — J). 

The remaining component of the no-jump dynamics is 
the driving Hamiltonian. This is the term that ensures 
that only the population in the steady state subspace will 
increase. Indeed, it affects all states with J ^ 0, trans- 
ferring population from the J z = — J 7^ into states with 
norm decay, while leaving the subspace of steady states 
unaffected. Therefore, the net effect of the dynamics be- 
tween detection events is to increase the population in 
the dark subspace. 

However, this dynamics is not sufficient to protect a 
particular highly entangled steady state of the system. 
As described in Section |III A[ the dark subspace of the 
dynamics has multiple linearly independent states. These 
states evolve exactly in the same way under the no-jump 
evolution, which will not be able to select one of these 
states over the others, as it can for states with J ^ 0. 
The fact that the system can be in different dark states 
in the quantum trajectory picture means that in the full 
density matrix description the state is in a mixed state, 
which has lower entanglement than the pure dark states. 
To ensure that only one target state in the subspace of 
steady states is accessed, we need to make sure that the 



system can't reach the other orthogonal states in this 
subspace. 



V. GENERAL FEEDBACK STRATEGY 

A general strategy for achieving a pure steady state can 
now be formulated. We need to find a feedback unitary 
U that achieves the two tasks described above: i) elim- 
inating the unwanted states in the dark subspace, and 
ii) providing transitions between states so that there is a 
path to the desired state. This will guarantee that the 
system will dynamically evolve towards a single steady 
state from any non-steady initial state. 

To design a feedback fulfilling these requirements let 
us first chose a target state from the set of steady states 
| J = 0, J z = 0,A§). We will label this state by A{j. All 
the states orthogonal to the target state in this subspace 
will then form the set of unwanted dark states |^ u ) = 
| J = 0, J z = 0, Aq 7^ Aq). These unwanted states can 
only be accessed from other subspaces via the combined 
jump operator J7J_. A way to avoid these transitions is 
to have the feedback unitary U connecting to \i/j u ) from 
a state |^) that cannot be reached through J_, i.e. , 



(a) 



^3 



J- 



(12) 



for any state The most obvious states for this are the 
unwanted dark states themselves so that U\ij) n ) = 
although any state with J z = +J would do (as there are 
no states for which J-\ip) has J z = +J). As a result, all 
other matrix elements in U connecting to the unwanted 
states will be 0: 



(^ u \U\J u J z ^^X k j)=0 



(13) 



for all i. 

The second task is accomplished by engineering a U 
that mixes all the subspaces with J / 0. This unitary 
is such that the jump operator UJ_ induces connections 
between the different J subspaces and provides to each 
one of them a path to the target state. Figure [3] shows 
different schemes to establish these connections. Suppose 
that a state in one space and ^2) in another are 
such that (ipi\U\ip2) 7^ (double-arrowed connection in 
Fig. [3]). If l^i) has J z ^ + J then there exists a state ^3) 
in the same subspace such that J-\ips) — (dashed 
arrow), so that UJ-\ips) = ^2) (solid arrow). Under the 
action of jumps, the system can then transition from the 
first subspace to the second. If ^2) also has J z ^ -f J, 
then the transition can happen in both directions (as 
shown by the solid arrows in Fig. [3]-b), otherwise the 
transition is unidirectional as depicted in Fig. [3}a. The 
path towards the target state can then be built by a chain 
of similar connections between subspaces. 

The general strategy presented in this section has a lot 
of freedom, with a variety of possible feedback operators 
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FIG. 3: Transitions between subspaces (columns) induced by 
the feedback operation. The grey doublepointed arrow shows 
the interaction via the unitary U, the dotted arrow shows the 
effect of the J_ operator, while the solid black arrow shows 
the net effect of the combined UJ- operation. Diagrams show 
transformations that results in jumps that can go from the left 
to the right subspace (a) (one-way transitions), or in both 
ways (b) (two-way transitions). 
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FIG. 4: Diagram of the angular momentum basis states for 
a four-partite system. States in the same J subspace are 
grouped together with the vertical ordering following the dif- 
ferent values of J z . Different columns within the same J sub- 
space are labelled with different Aj. 



fulfilling the requirements for stabilisation of a single en- 
tangled steady state. These choices may differ, however, 
on many other aspects such as the time it takes to reach 
the steady state, the robustness against imperfections or 
natural decoherence processes, and also on how simple 
it would be to experimentally implement them. We will 
leave these issues to the next section where we discuss the 
concrete case of controlling a system with four particles. 



VI. EXAMPLE: FOUR-PARTITE SYSTEM 

The subspace structure for a 4-partite system is shown 
in Fig. [4j States (represented by circles) within the same 
J subspace are grouped together, with the vertical posi- 
tion indicating the different J z projections. Although the 
configuration in this picture is general, there is some free- 
dom in the basis such that the choice of specific states 
associated to each circle is arbitrary. If we choose two 
basis states and ^2) for the J = subspace, for 
example, then any linear combination of these states will 
give a different valid choice of basis. 

This schematic representation of the subspaces is very 
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FIG. 5: Schematics of two different feedback schemes. Di- 
agrams (a) and (b) show the connections between states in- 
duced by U while (c) and (d) represent the effect of UJ-. The 
coupling given by the driving Hamiltonian is schematically 
shown in (e). One scheme (a and c) has one-way transitions 
between J subspaces, and the other (b and d) has two-way 
transitions that still provide a path to the target state. 



useful to analyse the effect of feedback. Two choices of U 
that follow the strategy described in Section [V] are shown 
in Figure pi (a and b), along with a representation of U J_ 
for each (c and d). Connecting arrows represent pairs of 
states that are coupled through the operators shown in 
the figure. One feedback uses two-way transitions that 
don't involve states with J z = +J (b and d in Fig. [5|, 
and the other involves one-way connections in a path 
from the space with J = 2 through to the steady state 
(a and c in Fig. [5J. Within each invariant irreducible 
subspace {J, Xj} the different spin projection states are 
mixed through the Hamiltonian term, and this is shown 
in Fig. [5}e. 

In both schemes we can see the requirements described 
in Section |V| fulfilled: The feedback unitary connecting 
the unwanted dark state ^ u to itself such that it can't be 
populated from other states, and a path connecting the 
different subspaces to the target state ijj t . In the direct 
one-way scheme of Fig. [5}c the path towards the target 
state is obvious and the Hamiltonian interaction depicted 
in Fig. [5}e ensures that the system will not be trapped 
in any of the J z = — J states (lowest circles in the J = 2 
and J = 1 subspaces). 

Despite leading to exactly the same final target state, 
these two choices of feedback are radically different dy- 
namically. Feedback schemes with unidirectional connec- 
tions are more efficient, in the sense that they bring the 
system to the target state faster. Indeed, as the system 
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FIG. 6: Simulations showing the effectiveness of the two feed- 
backs shown in Figure [5] with a spontaneous emission rate of 
10 _3 r. Each curve shows the overlap of the system state 
with a set of basis states. The thin lines represent the result 
of a single trajectory simulation and show the overlap with 
J = 2 (solid grey line), J = 1 (dotted blue line) and the tar- 
get steady state (solid red line). The thick red line shows the 
overlap with the target steady state for a full master equation 
simulations. The one-way feedback (a) takes fewer jumps and 
shorter time to reach the steady state than the two-way feed- 
back (b). The dips indicated by the arrows in each of these 
simulations represent spontaneous emission events. 



moves from the J = 2 to the J = 1 subspace via a jump 
process, it can't go back. From the J = 1, a detection 
event followed by the feedback application can only move 
the system within the subspace or to the dark (J = 0) 
subspace. However, in the case of bidirectional connec- 
tions, once in the J = 1 subspace the system can either 
follow the path towards the J = or go back to the J = 2 
subspace. This difference in efficiency can be seen in Fig- 
ure [6j where the overlap between the state of the system 
and the target state as a function of time is plotted for 
a single quantum trajectory. In the unidirectional case 
(Fig. |6]-a) the system reaches the steady state, on aver- 
age, ten times faster and with ten times less jumps than 
the two-way feedback showed in Fig. [6}b. 

The issue of time efficiency turns out to be of cru- 
cial importance when a natural decoherence effect, such 
as spontaneous emission, is taken into account. Sponta- 
neous emission is described by the third term in Eq. (3 
and a single decay event (indicated by arrows in Fig. 6 
will knock the system out of the dark subspace. After 
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the emission, the feedback control dynamics will have to 
restart the process of pushing the system to the target 
state. The effectiveness of the control then becomes a 
matter of competing feedback and spontaneous emission 
rates. In the simulations shown in Fig. |6j this means 
that the unidirectional feedback will be not only faster 
than the two-way feedback to recover after a spontaneous 
emission event, but also more effective as the systems will 
spend more time in the target steady state. The conse- 
quence of that for the unconditioned evolution, i.e. the 
solution of the master equation for the density matrix of 
the system, is that the entanglement of the system with 
spontaneous emission will be closer to the ideal value for 
the unidirectional feedback. This behaviour can be seen 
from the overlap of the unconditioned state p(t) with the 
target state shown by the smooth bold lines in Fig. [6j 



A. Simple implementations of the feedback 

In the examples presented above, we discussed how to 
construct feedback unitaries that satisfy the conditions 
for state preparation and protection, their efficiency and 
efficacy against spontaneous emission. So far we haven't 
considered the specific form of these feedback operations 
in terms of physical processes happening to the atoms. In 
principle, a generic U obeying our feedback strategy can 
be quite complicated, even involving multiple qubit inter- 
actions. Given that in real experiments one would have 
limitations to produce a feedback unitary U, in this Sec- 
tion we will change the focus and ask a different question: 
can we find feedback forms that satisfy (maybe only ap- 
proximately) our general requirements and that are yet 
simple to implement? We will answer this question by 
providing an example of a feedback that is based only on 
single atom operations and that is able to approximately 
protect one type of entangled state in the dark subspace 
of the problem. 

The steady state subspace of the four-partite system 
includes any state that is the tensor product of two anti- 
symetric Bell states as, for example, 



1 



\4>bb) = ^ (\ge)i2 



- \eg)i2) 
\geeg) - 



1 



\egge) 



(\ge)s4 - \eg)s4) 
+ \egeg)). (14) 



This state corresponds to the choice a = = <fi = 
and /3 = 1/2 in Eq. (|6|, and this is the state we can 
eliminate from the dynamics. As these antisymmetric 
Bell states are eigenstates of the two-partite total spin 
with J 12 = J34 = 0, they are dark to any local operators 
that are symmetric over the exchange of the two particles 
in each pair. Let's consider a feedback Hamiltonian of the 
form 



which corresponds to driving each atom with different 
intensities. If the coefficients ai are chosen such that a\ — 
d2 and a 3 = «4, then F\^bb) = and, consequently, 

we have U\^bb) = ^ %f \^bb) = \^bb), as required for 
the elimination of unwanted steady states described in 
Section [VJ 

The target state will now be the state with J = that 
is orthogonal to \^bb)- This state in this case is 



(2\ggee) + 2\eegg) - \gege) 
-\geeg) - \egge) - \egeg) ) , 



(16) 



which is a highly entangled 4-qubit singlet [38] . 

Figure [7] shows the overlap of the steady state with 
the target state \ip t ) for the choice a\ = CL2 = A and 
as = CL4 = —A, for various values of A. Without spon- 
taneous emission (dashed line) the system goes to the 
steady state (excluding the values A = kn with k in- 
teger). However, when spontaneous emission is taken 
into account the overlap is less than 80%. The reason 
for that is that our choice of a^s only achieves the first 
part of the strategy: eliminating all but one state in the 
dark subspace. A feedback of this form does not achieve 
the second goal of mixing the other states, as it leaves 
some states in the J = 1 subspace closed under the dy- 
namics. Without spontaneous emission the mixing of all 
subspaces is less important since the system will only 
occupy an isolated subspace if the initial state of the sys- 
tem (in our case \gggg)) has a component in this space. 
Spontaneous emission, however, promotes transitions to 
states that are not coupled via feedback and the system 
can be trapped in those states, hence not reaching the 
final target. Note that with spontaneous emission the 
initial condition is not important as there can be only 
one steady state [32] . 

To circumvent this problem and allow a mixing of all 
the other subspaces, we introduce a small deviation from 
the conditions a\ = (12 and a% = 0,4 by considering the 
following feedback Hamiltonian: 

F = Ao\ + A{1 - e)5f - Aa% - A(l - e)&%. (17) 

Figure [8] shows the overlap between the final steady state 
and the target state as a function of A and e with the 
spontaneous emission rates 7i = 7 = 10 _3 r. We see that 
the introduction of the parameter e is enough to connect 
the subspaces that were previously isolated mitigating, in 
this way, the effects of spontaneous emission. For e 7^ 
and A ^ kir the overlap is close to unity, demonstrating 
the efficiency of this particular local feedback to generate 
the many-body singlet of Eq. (16). 



i=l 



(15) 



VII. CONCLUSIONS 

In this paper we have shown how pure entangled states 
of a system of multiple atoms in a cavity can be prepared 
and protected using quantum-jump-based feedback. We 
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FIG. 7: Overlap between the fi nal s tate and the target steady 
state for the feedback of Eq. (15) with a\ — a<i — A and 
C13 = a a — — A as a function of A. The simulations show 
results without spontaneous emission (dotted blue line) and 
with 7 = T x 10 -3 (solid red line). For both of these simula- 
tions the system started with all atoms in the ground state. 
Spontaneous emission induces transitions to trapped states 
that do not connect to the target state, considerably decreas- 
ing the overlap. 




Overlap 



FIG. 8: Overlap between the final state and the target steady 
state for the feedback in Eq. (17) as a function of A and 
e. The spontaneous emission rates for these simulations is 
7 = r x 10 -3 . As in Fig. [7] there are dips in the overlap when 
A = kir. For e / the overlap gets closer to one as compared 
to the curve in Fig. [7] 



have formulated a general strategy that explores the 
steady state properties of the master equation and how 
they can be modified via feedback control. This engi- 
neered evolution dynamically steers the system towards 
a unique target steady state, inducing, in this way, a ro- 
bust generation of multipartite entangled states. 



The general requirements for an effective feedback pre- 
sented here are independent of the number of atoms in 
the system. For this reason, the feedback operators sat- 
isfying these conditions can, in principle, be constructed 
for any number of qubits. However, practical consid- 
erations would impose extra constraints on the form of 
the unitary U that could be reasonably implemented in 
real experiments. Here we partially addressed this is- 
sue by showing that a simple feedback based on single- 
atom operations can be used to stabilise a given class of 
multipartite entangled states in the case of four qubits. 
This example, together with the freedom in the choice of 
the feedback J7, indicate that it may be possible to find 
equally simple feedback Hamiltonians for systems with 
arbitrary number of constituents. At the moment, a gen- 
eral approach to this problem remains an open question. 



Finally, we should mention that the dark subspace can- 
not be changed by the feedback as it is solely determined 
by the measurement operator. Any attempt to modify 
the range of accessible pure steady states would then 
require engineering a new measurement scheme on the 
system. The ability to do this, however, depends on the 
specific details of the model and is subject to the limi- 
tations imposed by the physical system analysed. Here 
we have shown that, for a specific measurement scheme, 
quantum-jump-based feedback can be used to alter the 
dynamics of the system and the stability of steady states. 
In the future, it would be interesting to consider new 
measurement procedures that, allied to feedback design, 
would open up the possibilities to engineer a broader 
range of interesting multipartite steady states. 
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